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Abstract 

Three-loop corrections to the scalar and pseudo-scalar current correlator are calcu- 
lated. By applying the large momentum expansion mass terms up to order (m 2 /g 2 ) 4 
are evaluated analytically. As an application 0(cy 2 ) corrections to the decay of a 
scalar and pseudo-scalar Higgs boson into top quarks are considered. It is shown 
that for a Higgs mass not far above the ti threshold these higher order mass correc- 
tions are necessary to get reliable results. 
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1 Introduction and notation 



A crucial question in elementary particle physics is whether nature makes use of the 
spontaneous symmetry breaking for generating the masses or not. In the minimal standard 
model (SM) this mechanism requires the existence of a scalar particle, the Higgs boson. 
Extensions of the SM, e.g. models with more than one Higgs doublet or supersymmetric 
versions of the SM, predict also pseudo-scalar Higgs bosons, A. 

However, up to now there is no experimental evidence for such a particle and the direct 
search at LEP via the process e + e~ — > ffH rules out the mass range Mh < 65.6 GeV 
with a 95% confidence level (CL) [|IJ. Assuming the validity of the SM the precision data 
are sensitive to the Higgs boson and a recent global fit yields Mh = 149tg| 8 GeV together 
with a 95% CL upper bound of 550 GeV (in this context see also 0). Theoretical 
arguments based on unitarity or the validity of perturbation theory request an upper limit 
on the Higgs mass of about 1 TeV. 

In this paper we consider the decay of a Higgs boson, scalar or pseudo-scalar, with 
mass above 400 GeV, into a top-antitop quark pair. Since the velocity of one of the quarks 
for this mass range is v J> 0.5, threshold effects can safely be neglected and the decay rate 
can be treated in a purely perturbative way. The full mass dependence for this process is 
only known to 0(a s ) while to order a 2 only mass terms up to M 2 /M H , A are available 
|5|, H [7|, |S|, where M t is the top quark mass, and they show up to be quite sizeable |J. 
0(af) corrections have recently been evaluated in the massless case 0. 

In the following, making use of a technique recently developed for the automatic com- 
putation of mass corrections to the vector current correlator [[10], |i~lf , we will demonstrate 
that the inclusion of mass terms up to order (M 2 /M H , A ) 4 leads to reliable results at 
0(a 2 ) in the considered mass range. We note that the results presented in this paper for 
the decay into top quarks may be generalized to any fermion species. However, even for 
bottom quarks already the quadratic mass corrections are very small. 

To fix the notation we define: 



where n s (g 2 ) (n p (g 2 )) represents the scalar (pseudo-scalar) current correlator in momen- 
tum space. The currents are given by 



It should be noted that as far as renormalization is concerned it is convenient to consider 
the combination mj s in the on-shell scheme and mj s in the MS scheme, where m (m) is 
the generic pole (MS) mass, in order to avoid additional renormalization constants. The 
physical observable R 5 (s) is related to U s (q 2 ) via the relation 
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The current correlator can be written as 
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and similarly for R 5 (s). The normalization in Eq. (j3|) guarantees that R^ ,s (s) — > 3 
for s — > oo. The colour factors (C^ = (N% — 1)/(2N C ) and Ca = N c ) correspond to 
the Casimir operators of the fundamental and adjoint representations, respectively. For 
the numerical evaluation we set N c = 3. The trace normalization of the fundamental 
representation is T = 1/2. The number of light (massless) quark flavours is denoted by 

m. 

In Eq. ([]) n^' 5 is the abelian contribution which also exists in QED, and IT^ 5 results 
from the non-abelian structure specific for QCD. The contribution of diagrams containing 
a second massless or massive quark loop is denoted by Elj and Up^' S , respectively. 
represents the terms arising from the double-triangle diagram and is called the singlet 
contribution. The case for the double-bubble diagram where the inner quark mass is 
much heavier than the outer one is not listed as it does not contribute in cases of physical 
interest . 

The outline of the paper is as follows: In Section |2] the method is briefly described and 
the results for the correlator functions are given. In Section |3] mass effects for the decay 
of a scalar and pseudo-scalar Higgs particle into top quarks are presented and the results 
are discussed. 



2 Results for the current correlators 

In this section we keep the discussion general and consider the scalar and pseudo-scalar 
current correlators with a generic quark mass m. In Section ^| the specification to the case 
m = M t is performed and numerical values are given. 

Let us in a first step outline the main ideas of the large momentum expansion and 
its realization by a computer^. The large momentum procedure [12] requires the iden- 



tification of certain subgraphs associated with the Feynman diagram to be calculated. 
These subgraphs then have to be expanded in their small dimensional quantities (i.e. all 
except the large momenta) and the resulting terms, being products of tadpoles and mass- 
less integrals, have to be calculated. The number of these terms increases rapidly with 
the number of loops. In our case of a massive current correlator, the one- and two-loop 
diagrams generate 17 terms altogether, so that they still can be treated by hand. In the 
three-loop case, however, there are 19 topologies, contributing 266 terms, and their gen- 
eration by hand is not feasible. Therefore we completely automated the large momentum 
procedure for massive two-point functions and directly fed the generated terms to the 
FORM pi packages MATAD and MINCER 0]. After calculating the single terms the 



results of all diagrams were added giving a finite expression after renormalization. 

We are now prepared to present the result. Let us first consider the scalar correlator. 
The results for the different contributions read in the MS scheme (l qm = ln(— g 2 /m 2 ), l qfl = 
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where ( is Riemann's zeta function with values (2 = n 2 /G, (3 w 1.20206 £4 = 7r 4 /90 and 
C5 ~ 1.03693. The constant .84 is typical for three- loop tadpole integrals and is given by 
i?4 « —1.76280 [|lf|. The corresponding results in the on-shell scheme are obtained by 
substituting in the combination m 2 n(g 2 ) the MS mass, fh, w.r.t. the on-shell mass [1(J . 



For the pseudo-scalar current correlator a minor complication arises in connection with 
75. While the non-singlet contribution allows the use of the anticommuting definition, for 
the singlet diagram we take the prescription for 75 introduced in ||17|| . Thereby we follow 
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the strategy outlined in |L8], [7| and replace 75 in the vertices according to 

75 " 
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where 7^ p<T l represents the antisymmetrized product of four 7-matrices and may be writ- 
ten as 
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Because the double-triangle diagram contains no sub-divergencies and consequently has 
a finite imaginary part it is allowed to contract the new polarization tensor with eight 
indices with the product of four metric tensors and compute the scalar integrals in the 
same way as for the non-singlet contributions. Even more, in contrast to the scalar case 
U^' p is finite from the very beginning. 

The results for the pseudo-scalar case read: 

2 + 4 Iqm 



n(o),p 



16tt 2 



m I m 

4: — 2lqn + 4:—l q n + I — 



777 

Q 



2\ 3 



+ 1^ 



20 

g "I - 8 



/ 777 

+ k 



59 

2 1 20 ig m 



+ 



(15) 



n(i),p 



16tt 2 l 8 



+ 



711 



131 ^ 

2\ 2 r 



— o 
777 



17 3 2 



-6 + 12C 3 + 4L M -6/ 



q/l 



/ 777 

+ br 



2\ 3 



4 + 27 l qm + 6 /^ + (15 - 18 V) V 
515 628 82 2 

g + g •'qm + 2 ' 9 m + (52 48 L qm ) l qfM 



+ 



' - 2\ 4 r 
777 \ 



36985 3101 , 539 , 2 

1 / -I / -I- 

144 + 18 9m+ 6 qm + 



(f -150 V,)'. 



+ 



(16) 



n 



(2),P 



3 r 1613 
16tt 2 1 64 
777 2 T1697 



9 / 691 9 



105 It- ^ /? 



9 



/ 777 
+ '^2- 



48 

2\ 2 



129 



+ 54C 3 -6C 4 -70C 5 



885 



13 



27 



--24Cs)V- T C + 6 ^ 



- — + 104 Cs-24 C4 + 15C5 + 4S4 



/1253 



V 16 



33CaM 



h qm 



-I 7/ + 6/ d 



195 423 



/ —27/1/ + 



gm / "'QA' 



^ + ^/ )/ 2 



+ 



' - 2\ 3 
777 \ 



421657 
1944 



135C3-48C4 + 80C5 + 8S4 + 



/4481 



+ 52C 3 H. 



gm 



6 



839 l9 1223 ,„ 
-I 1 -\ Z - 

1 g "gm 1 27 <? m 

+ (-192 + 144 l qm )l 2 m 



/909 1028 
V~2 3~~ 



/ _1R4/ 2 )/ 



/m 2 \ 4 



™_| C3 _ 120C4 + 2 80 C5 + 20 B4 



+ 



+ 



/ 751555 

V 1296 
/70621 

V 32 ~ 



203 CaW 



b qm 



534067 2 
+ _ 864 _ 9m 



38753 



3124 _ 2695 2 . ; 



216 ^ 

6225 1125 
— ^- + 



2 l qm J l qll 



+ ... 



(17) 



n 



NA 



3 f 14419 75 9 5 
16^2 I 288 ~~ T C3 ~ 8 C4 ~~ 2 Cs + 



893 31 . , 



11/2 _11,3 

12 ^ 24 9/i 



— 2 

m 



+ f C3 + 3C 4 + 5C 5+ (^-17 C 3)^ 



_ 11^/2 ,11/3 

12 "* 6 w 

2\ 2 r 1417 



+ 1? 



96 



^ > > 5 A „ Z3925 
+ T C3 + 12C4 + -C 5 -25 4+ (^ 



38 2 11 3 

+ 2 \m + g 'gm + 
2\ 3 



49 
4 

397 



_ 49 / -H/ 2 ) 7 + ( - 

24 9,71 9 y" 1 / 9M 



+ 7" 



134501 63 



972 
505 



- y C 3 + 24C4 



2 

20 C 5 



55 



15 \ 

4 l ?m J V 



45 4 + 



/ 



497 ,2 

-I / 2 _i_ 

1 g qm 1 54 "gra 



/13231 

+ hor 



3473 451 2 

" ^gm — » 



27 



gm J 



( 170609 
V 648 

In 



-40C 3 W 



/m 2 \ 4 

+ 1? 



3008885 
3888 



823 
~24 



C3 + 60C4-70C5- 105 4 



/1489679 1/1CA \, 
+ (^592-- 145C3 ) Z - + 



74807 



+ 



288 m 

820055 77761 , 5929 , 2 
/ Z 

l gm 



8941 7 o 
+ — Z 3 



1728 



216 



Y2 i m J 



216 qm 
3905 



^ l 9» I 'g/i 



+ . 



(18) 



n 



16tt 2 

m 2 
+ — 



{-§ + 4C3+ (f - 2 ^)'»--y , «- + g4 



36 
31 



31 4 / 9 \ 7 2 2 3 
~36~3 C3+ r2 +4C3 J ^ + 3^~3^< 



7 



3 9/1 



7 



17 



53 257 10 2 2 3 



g qm 



+ (-y + 14/ gm + 2/ 2 
,2 



+ -IT 



12551 
486 



16 4970 

y C3 ~^i 9m 



g ~<jrm 
5 ■>, 

133 



Z 



-Z 2 -^Z 3 

g gm 2^ 9 m 



+ 



905 988 82 2 



,26 



qm I ^ 



(m 2 \ 4 



573913 4 13547 1315 2 287 3 

~ 3 81 m 24~~ m 54~ m 



3888 



+ 



58285 5351 



—I 

16tt 2 1 

— 2 

m 

+ L 36 
( m 2 \ 2 



432 
511 



+ 



54 



-Z 



539 2 \ 
+ 18 9m J M 



355 
T2~ 



25 Z„ m ) Z 



qm ) l qtl 



+4c 3 +(f -2c 3 )v-^/;„ + -/, 3 



8 



11 



36 



6 9/1 6 gfM 



701 293 

^ + 6C 3 -— U 



11 



3 "<2M 3 9/* 



-/ 2 --/ 3 

( — + 14 Zgm + 2 Z ?m j l qix + ( — — 3 Zgm, J Z g ^ 



/ m 



+ 



6 

2\ 3 



43319 596 3064 623 2 92 3 

1 rT~ ^ 3 n^v~ i m 27~ 9m — 27 m 



905 



486 
98£ 

+ 



9 



27 



27 27 



?™ j ^qn 



+ 



< - 2\ 4 



3892289 2333 . 20533 



7776 



+ 



9 



+ 



58285 5351 



Ca- 
539 



24871 



zL- 



545 



+ . 



z; 



72 qm 216 " m 27 gm 



432 



54 



-^g qm J <?M 



355 
~Y2 



25 Z om ) Z 



qm ) l qtl 



+ . 



167T 2 

+ 



2 r 



+ 
+ 



m 



- 2\ 2 



2\ 3 



,7 



2\4 



16C 3 -20Ce 



- 44 + 24C 3 + (-12-72C 3 )Z (?m 



f + H4 Ca + (~ ~ S6 Ca) km ~ 63 l\ m - 12 Z 3 



'qm 



68146 1288 



243 



+ 



9 



■Ca + 



( 



7727 



V 81 



-80CaK m -86Z 2 m -44Z 3 m 



8 



The results for the one- and two-loop case have been checked against the analytical exact 
expressions for the polarization functions [19] in the limit of large external momentum q. 



3 The decays H — ► tt and A 



tt 



As an application of the results of the previous section we consider the decay of a scalar 
and pseudo-scalar Higgs boson into a top-quark pair. Thereby we neglect all light quark 
masses and denote the MS top quark mass by m t and the on-shell one by M t . Let us 
in a first step review the Born result and the first order QCD corrections. Taking the 
imaginary part of Eqs. (§) and (0) and transforming the combination mfR(s) into the 
on-shell scheme (16[ one arrives at (L ms = ln(M 2 /s)) 
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for the scalar case where /i 2 = s is chosen. For the pseudo-scalar current we get from 
Eqs. (HD and 
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In Fig. |I] the results are plotted together with the exact expressions. For the abscissa 
the variable x = IMtjyfs is chosen. Already the approximations including only the 
quadratic terms provide a good agreement up to x ~ 0.4. It is, however, dramatically 
improved by incorporating higher terms in M 2 /s, leading to an excellent approximation 
almost up to x = 1 for R^°'' S ,R^' S and R^' p and up to x ~ 0.85 for This serves 

as a strong motivation to evaluate higher order mass terms for the 0{a 2 s ) corrections. 

At 0(a 2 s ) the results for R^' s (s), separated into the contributions from the different 
colour factors, read in the on-shell scheme: 
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Figure 1: R^' 5 and R^' s , 5 = s,p, plotted against x = 2M t /y/s. Successively higher 
order terms in (M 2 /s) n : Wide dots: n = 0; dashed: n — 1,2, 3; solid: n = 4; narrow dots: 
exact. 
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In the pseudo-scalar case we get: 
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The decay width of a scalar and pseudo-scalar Higgs boson is then given by 
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where -R^' 15 is the contribution from the pure gluonic cut appearing in the imaginary part 
of the singlet diagrams which has to be subtracted if we are interested only in the fermionic 
final states. Rfg' 5 is known analytically both for the scalar [2(| and pseudo-scalar case 
2j| . For completeness we list the expansion in M^/s: 



3{f 



2 + 6C2 + y C4 + (-1 + ~C 2 ) L 2 ms + ~Li 



-24C 2 - 90C 4 + (-4 + 6C2) L ms + (4 - 12Ca) L 



2 

ms 



3 i i 

m 

2\3 



' ms ms 



-4 + 6C2 + I8OC4 + (10 - 39C 2 ) L ms + (3 + 24( 2 ) L 



2 

ms 



13 /-3 ,or4 



2\ 4 r 



+ - 



V s 



Mf 



44 



22 



4 - 30(2 + ( y + 44C 2 J L ms - 15LL + y^L 



+ 



45 ^ + 3 a 2 +"y 



1 

ms 



+ 



' M 2 



6C2-^ms + ^L) 



(37) 



13 





(M 2 
\ « 


+ 




(M^ 
\ s 


+ 



6( 2 + 9( 2 L ms + 3L 2 ms + h 3 ms 
I8C2 + (4 + 20C 2 ) L ms + 9L 2 rns + ™L 3 ms 



+ ... . (3* 



Note that the analytical structure is not changed if M t is replaced by the MS mass because 
Rfg ,& only contributes at order a 2 . For the numerical results presented below we only use 
the inclusive quantity R s (s). 

The results for the constant and first mass corrections are in agreement with the 
literature both for the scalar || [7| and pseudo-scalar case || . A strong check is provided 
for the terms proportional to n\ where a successful comparison with the exact results for 
the scalar [^, |23|] and pseudo-scalar [^| contributions was possible^. 



However, one objection is in order here, namely the fact that R^' S , Rpf , R^' 5 and 



In Figs. and |3| the contributions of the different colour factors are plotted including 
successively higher orders in M 2 /s. The notation is the same as in Fig. |l| Motivated by 
the behavior of the Born- and one-loop terms of Fig. [I], one may consider as a measure 
of validity for each curve the range in x where it coincides with the one containing most 
powers of x (i.e., the solid line in Figs. and |3[). This criterion is further justified by 

(2) 8 

the light fermion contribution R\ , where also the exact result is known. For the scalar 
case, i?^' 5 , the agreement of the exact result and the M t 8 /s 4 -terms is almost perfect up to 
x — 1, while for R^ the convergence above x = 0.9 gets poorer, but the above argument 
concerning the validity holds perfectly in both cases. 
However, one objection 

(2) 8 (2) 8 

R s , in contrast to R\ , exhibit a four-particle threshold at x — 1/2 which may spoil 
convergence for x > 1/2. Indeed, the curves for R^p' 5 , R^' S (Fig. ^) develop a relatively 
large spread in this x-range. Nevertheless, the range of validity as defined above extends 

( 2) 8 

to x = 0.7 — 0.8. Despite this objection, the curves for R Nj ( seem to converge very well 

(2) 8 

again almost up to x = 1, those for R A at least up to x ~ 0.85. In particular for the 
scalar case, where the exact result is known to be zercQ at x — 1, the behavior of the 
approximations is quite promising. 

Summarizing, one can see that while the quadratic mass terms seem to reasonably 
approximate the result up to x = 0.3 — 0.4 for the scalar and x ~ 0.5 for the pseudo-scalar 
case, the quartic terms considerably improve the expansions in most cases up to values 
of x ~ 0.85, in some cases even up to x = 1. The influence of the higher order terms is 
rather small and practically only needed for R^' 5 and Rg' S . 

It is interesting to compare the results both for the on-shell and the MS-definition of 
the quark mass. Therefore we set /x 2 = s which is the natural choice because it eliminates 
all logarithms of the constant and quadratic terms in the MS scheme. Starting from the 
quartic term, however, In M 2 / s terms remain, resulting from non-trivial operators which 



are absent before 24 . In contrast to the MS scheme the results in the on-shell scheme 



develop In M 2 / s terms starting from the lowest order. It is nevertheless instructive to 



2 We would like to thank the authors of |23) for providing us with the results prior to publication. 

.(2; 
s 



3 R < g ^ because the gluon cut is still present. 
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Figure 2: R^' S with 5 = s,p and n = A, NA, I plotted against x. The same notation as 
in Fig. ^ is adopted. 



compare the expansions both in Mf / s and a s (s) in the two schemes. Therefore the 
renormalization group invariant quantities m 2 (s) R s (s) / M 2 and R 5 (s) are considered. We 
choose rii = 5, af\hl 2 z ) = 0.118, M t = 175 GeV and M H/A = 450 GeV which corresponds 
to x ~ 0.78. The corresponding MS top mass, m t (M H /A), evaluates to 77^ (450 GeV) = 
155 GeV and the two-loop beta function leads to a s ((450 GeV) 2 ) = a( 6 )((450 GeV) 2 ) w 
0.096. In Tab. [1] the numbers are listed. It can be seen that the quartic terms are still 
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x x 

Figure 3: R^' s with 5 = s,p and n = F, S plotted against x. The same notation as in 
Fig. D is adopted. 



sizeable and in many cases numerically of the same order as the constant and quadratic 
ones. They become particularly important when the first two terms happen to cancel 
to a large extent. The higher order terms give only small contributions at least in the 
three-loop case. 

Usually in the MS scheme a faster convergence of the QCD series is expected in the 
sense that the coefficients in front of (a s /7r) n are smaller than in the on-shell scheme. 
Considering the sum of all approximation terms available this happens to be true only 
for the scalar correlator whereas in the pseudo-scalar case it is vice versa. A possible 
explanation might be that the coefficient of the leading term in the on-shell scheme is 
already rather small because of an accidental cancellation between the constant and the 
logarithm. For the scalar correlator there is in addition a large cancellation between the 
first three terms in the MS scheme so that the sum turns out to be also rather small. 
This is not the case for RP where the first three terms have the same sign. However, the 
scheme dependence is drastically reduced by the inclusion of the 0{a 2 s ) terms as can be 
seen in Tab. [I]. 

Finally we want to present numerical values for the physical decay rate into top quarks. 
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(M t 2 )° 


(M?) 1 


(M t 2 ) 2 




(M t 2 ) 4 
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exact 




scalar, on-shell 








1.000 


-0.907 


0.137 


0.014 


0.003 


0.247 


0.248 


C F R {1) < s /3 


-0.778 


5.646 


-3.080 


0.010 


0.004 


1.802 


1.802 


RW' s /3 


-35.803 


77.056 


-17.792 


-5.347 


-0.680 


17.435 


— 




0.943 


-0.663 


0.026 


0.009 


0.003 


0.318 






scalar, MS 






|^i? (0) <73 
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0.784 


-0.558 


0.066 


0.005 


0.001 


0.298 


0.299 




4.445 


-3.723 


-0.238 


0.143 


0.032 


0.659 


0.673 


1 

^i? (2) 'V3 


21.799 


-7.606 


-15.334 


0.680 


0.546 


0.086 


- 


S,(a s /7r) 1 


0.941 


-0.679 


0.045 


0.010 


0.002 


0.319 






pseudo-scalar, on-shell 






fl(o),P/ 3 


1.000 


-0.302 


-0.046 


-0.014 


-0.005 


0.633 


0.629 




-0.778 


3.495 


0.417 


0.047 


0.027 


3.208 


3.238 




-35.803 


25.024 


12.173 


3.780 


1.293 


6.467 






0.943 


-0.172 


-0.022 


-0.009 


-0.003 


0.737 






pseudo-scalar, MS 








0.784 


-0.186 


-0.022 


-0.005 


-0.002 


0.569 


0.569 




4.445 


-0.248 


-0.215 


-0.119 


-0.043 


3.821 


3.791 


I"i R(2),p/3 


21.799 


9.854 


2.455 


-0.781 


-0.520 


32.807 




S l (a s / 7 r) 1 


0.941 


-0.185 


-0.026 


-0.010 


-0.003 


0.717 





Table 1: Numerical results for R s and RP both in the on-shell and MS scheme. The 
contributions from the mass terms (Mf) 1 , their sum (E) and, where available, the exact 
results are shown. Ej(a; s /7r) 1 is the sum of the 1-, 2- and 3-loop terms. The numbers 
correspond to M t = 175 GeV and M h /a = 450 GeV. The renormalization scale /x 2 is set 
to 8 = M* /A . 

For convenience the square of the on-shell mass, M 2 , is factored out both in the on-shell 
and MS scheme. Using Eq. ([36|) the expansions of the decay width for the scalar and 
pseudo-scalar case look like: 

GfMhM? 

r(H^tt) = jL 1 [0.7404 + 0.1652 + 0.0469 + 0.0019 - 0.0022], (39) 

4y 2n 

GfMhM? 

T(H^tt) = " 1 [0.8954 + 0.0604 - 0.0011 + 0.0014 - 0.0020], (40) 

4\/27r 
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T(A^ti) = GpM u M t [ L898 2 + 0,2941 + 0.0146 + 0.0035 - 0.0051] , (41) 
4v27T 

T(A^ti) = —^^^[1.7084 + 0.3503 + 0.0903 + 0.0016- 0.0038], (42) 
4y27r 

where the first and second number correspond to the Born and 0(a s ) correction. The 
0(a 2 s ) terms are separated into the non-singlet (third number) and singlet (fourth number) 
contribution. The last number corresponds to the gluonic cut which has to be subtracted. 
The last two numbers cancel each other to a large extent which means that the two- 
gluon-cut dominates the imaginary part of the double-triangle diagram. Furthermore one 
observes that the two-loop QCD corrections amount up to ss 5%. 

To conclude, we have computed analytically the first five terms of the scalar and 
pseudo-scalar current correlator in the expansion for large external momentum q. As 
an application the decay of a scalar and pseudo-scalar Higgs boson into top quarks is 
considered. In this case the higher order mass corrections, especially the quartic terms, 
emerge to be important. 
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